Abstract—Intergenerie hybridization 
between the epinepheline serranids Ce- 
Phalopholis fulva and Paranthias fur- 
cifer in waters off Bermuda was inves- 
tigated by using morphological and 
molecular characters. Putative hybrids, 
as well as members of each presumed 
parent species, were analyzed for 44 
morphological characters and screened 
for genetic variation at 16 nuclear 
allozyme loci, two nuelear (mDNA loci, 
and three mitochondrial (mt)DNA gene 
regions. Four of 16 allozyme loci, cre- 
atine kinase (CA-B*), fumarase (FH"), 
isocitrate dehydrogenase (/CDH-S*), 
and laetate dehydrogenase (LDH-B*), 
were unique in C. fulva and P. furcifer. 
Restriction fragments of two nuclear 
DNA intron regions, an aetin gene 
intron and the second intron in the 
S7 ribosomal protein gene, also exhib- 
ited consistent differences between 
the two presumed parent species. Re- 
striction fragments of three mtDNA 
regions—ND4, ATPase 6, and 125/168 
ribosomal RNA—were analyzed to 
identify maternal parentage of puta- 
tive hybrids. Both morphological data 
and nuclear genetie data were found to 
be consistent with the hypothesis that 
the putative hybrids were the result of 
interbreeding between C. fulva and P 
furcifer. Mean values of 38 morphologi- 
ca] eharaeters were different between 
presumed parent species, and putative 
hybrids were intermediate to presumed 
parent species for 33 of these charae- 
ters. A principal eomponent analysis 
of the morphological and meristie data 
was also consistent with hybridiza- 
tion between C. fulva and P. furcifer. 
Thirteen of 15 putative hybrids were 
heterozygous at all diagnostic nuelear 
loci, consistent with F, hybrids. Two 
putative hybrids were identified as 
post-F, hybrids based on homozygosity 
at one nuelear loeus each. Mitochon- 
drial DNA analysis showed that the 
maternal parent of all putative hybrid 
individuals was C. fulva. A survey of 
nuelear and mitochondrial loci of 57 C. 
fulva and 37 P. furcifer from Bermuda 
revealed no evidence of introgression 
between the parent species mediated 
by hybridization. 
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Poey (1860, 1875) described the genus 
Menephorus for two species, M. dubius 
Poey, 1860 and Af. punctiferus Poey, 
1875, each based on a single specimen 
of grouper that appeared to be interme- 
diate in morphology between the coney, 
Cephalopholis fulva and the creole- 
fish, Paranthias furcifer. Smith (1966) 
analyzed Poeys specimens for 45 
meristic and morphometric characters 
and found the specimens to be interme- 
diate to P. furcifer and C. fulva for 40 
of these characters—results consistent 
with intergeneric hybridization. Smith 
(1966) also noted that P furcifer may 
be more closely related to the epineph- 
eline serranids (such as C. fulva) than 
previously thought. Interest in this 
possible case of intergeneric hybridiza- 
tion was renewed in 1993 when Bermu- 
dian ichthyology student James Parris 
Jr. asked one of the authors if he would 


be interested in a specimen his father 
caught that “had the head of a coney 
and the tail of a barber (creole-fish).” 
Since that time several other putative 
hybrids have been caught off Bermuda 
(Smith-Vaniz et al., 1999), and we have 
initiated a program to collect putative 
hybrids from local fishermen to further 
investigate this phenomenon. 

The presumed parents of the puta- 
tive hybrids are members of different 
serranid genera that have strikingly 
different morphologies, occupy different 
ecological niches, and have different be- 
haviors (Heemstra and Randall, 1993; 
Smith, 1966). The major distinguishing 
features of P. furcifer are its forked 


* Contribution 2469 of the Virginia Institute 
of Marine Science, College of Wiliam and 
Mary, Gloucester Point, VA 23062-1346. 
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Figure 1 


Photograph of the creole-fish (top), Paranthias furcifer, the coney (bottom), Cephalopholis fulva, and the puta- 
tive hybrid (middle), Bermuda. (Courtesy of Jan Cordes.) 


caudal fin and the size and shape of its mouth (Fig. 1). 
Paranthias furcifer is dark red in color; darker dorsally, 
lightening ventrally (Heemstra and Randall, 1993). There 
is an orange spot at the upper end of the base of the pecto- 
ral fin and three white spots dorsal to the lateral line. Ber- 
muda is the northern limit of the distribution of P furcifer 
in the western Atlantic and the species occurs throughout 
the Bahamas and Antilles and along the American coast 
from the Gulf of Mexico south to Brazil (Smith, 1971; 
Smith-Vaniz et al., 1999). 

In contrast, C. fulva is usually scarlet in color and cov- 
ered with light blue-green spots, each surrounded by a 
black ring (Heemstra and Randall, 1993). There are two 
black spots on the edge of the lower jaw as well as on top 
of the caudal peduncle. Cephalopholis fulva has rounded 
caudal, anal, and dorsal fins, similar to fins of other 
epinepheline serranids (Heemstra and Randall, 1993) 
(Fig. 1). In the northwestern Atlantic, C. fulva has a distri- 
bution similar to that of P. furcifer; it occurs as far north 
as Bermuda, throughout the Bahamas, Antilles, and along 
the east coast of the Americas from South Carolina to 





Brazil (Heemstra and Randall, 1993; Smith, 1971; Smith- 
Vaniz et al., 1999). 

Morphologically, the putative hybrids are almost exactly 
intermediate to the parent species (Fig. 1). In his 1966 re- 
view, Smith noted that the hybrids have some characters 
unique to C. fulva and P. furcifer. For example, the puta- 
tive hybrid individuals have both a moderately forked tail 
and blue spots surrounded by a black ring. Paranthias 
furcifer is the only Atlantic grouper with a forked tail sug- 
gesting that it is one of the putative parents and C. fulva is 
the only Atlantic grouper that has blue spots with a black 
ring. Presence of both traits together in a single individual 
strongly suggests interbreeding between the two species. 
The geographic extent of the putative hybrid is not well 
known; however specimens exist from Cuba and Bermuda 
(Smith-Vaniz et al., 1999), and it has been reported from 
Jamaica (Thompson and Munro, 1978). 

Hybridization has traditionally been detected by using 
morphological characters, but increasingly, genetic analy- 
ses have also been used for this purpose. Allozyme elec- 
trophoresis provides a rapid and cost-effective method to 
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assess hybridization (Campton, 1987). Through the analy- 
sis of multiple loci it is possible to identify F, and post-F, 
hybrids, as well as to detect the introgression of alleles 
between species by hybrid backcrossing. It became evident 
however, that analysis of nuclear DNA loci is necessary 
in studies of hybridization as a means to overcome some 
of the sampling restrictions of allozyme analysis and to 
provide a survey of biparentally inherited genes (Verspoor 
and Hammar, 1991). Mitochondrial DNA (mtDNA) has 
also been used extensively in studies of hybridization. Be- 
cause mtDNA is maternally inherited in fishes, analysis of 
the molecule allows one to identify the maternal parent of 
F, hybrids, as well as the sexual preferences of F, hybrids 
and their offspring (Dowling et al., 1996). 

The key to using molecular markers to identify hybrid- 
ization is to find multiple independent nuclear loci and 
a mitochondrial gene region that have unique alleles in 
each putative parent species (Dowling et al., 1996). An F, 
hybrid would be heterozygous at all nuclear loci and have 
a mitochondrial haplotype identical to one parent species 
(Campton, 1987; Dowling et al., 1989). A backcrossed indi- 
vidual would be heterozygous at some diagnostic nuclear 
loci and homozygous at others. Therefore the power of 
demonstrating an F, hybrid, as opposed to a backcross 
or pure parent individual, increases with the number of 
nuclear loci examined. 

In this study, genetic information from four diagnos- 
tic allozyme loci, two diagnostic nuclear DNA loci, and 
three diagnostic mtDNA gene regions, was used to assess 
hybridization and introgression between Cephalopholis 
fulva and Paranthias furcifer in Bermuda waters. 


Materials and methods 


A total of 51 Cephalopholis fulva (Linnaeus, 1758), three 
C. cruentata (Lacepède, 1802), and 37 Paranthias furci- 
fer (Valenciennes, 1828) were collected from Bermuda 
with baited handlines or rotenone solution. In addition, 
six C. fulva and two C. cruentata were sampled from 
Navassa Island. Fifteen putative hybrids were captured 
by Bermudian fishermen using handlines or lobster traps. 
Cephalopholis cruentata was in-cluded in the study as a 
possible parent species of the putative hybrid and three 
Epinephelus guttatus (Linnaeus, 1758) specimens from Lee 
Stocking Island, Bahamas, were used in the preliminary 
mitochondrial DNA study. 

Specimens were frozen upon capture and stored at — 
20°C or -80°C and transported to the laboratory for analy- 
sis. For mitochondrial and nuclear DNA analysis, muscle 
tissues were removed and placed in storage buffer (0.25M 
EDTA, 20% DMSO and saturated with NaCl). For allo- 
zyme analysis, 1.5-cm® pieces of liver and muscle tissue 
were separately homogenized in 250 uL of chilled (4°C) 
grinding buffer (0.1 M Tris, 0.9 mM EDTA, and 0.05 mM 
NADP*, pH 7.2). Samples were centrifuged for 3 min at 
16,000 xg and stored at -80°C or analyzed immediately. 

Genomic DNA was isolated from a 1.0-cm® piece of mus- 
cle tissue by using the phenol/chloroform protocol of Win- 
nepenninckx et al. (1993) with the following modifications. 


Table 1 


National Museum of Natural History (USNM) and Vir- 
ginia Institute of Marine Science (VIMS) catalogue num- 
bers for specimens used in the morphological] analysis. 


Standard 


Species Collection Specimen length (mm) 


USNM 88717 232 
USNM 53134 217 
USNM 53134 176 
USNM 133689 Ail 
VIMS 10413 222 
VIMS 10414 213 
VIMS 10415 210 
VIMS 10416 170 
VIMS 10417 146 
USNM 320539 240 


P. furcifer VIMS 10410 274 
VIMS 10411 259 
VIMS 10412 274 
USNM 107108 218 
USNM 358541 237 
USNM 33255 207 
USNM 12540 185 
USNM 65605-8244 201 
USNM 65605-8246 176 
USNM 65605-8246 129 


Putative hybrid VIMS 10403 190 
VIMS 10402 199 
VIMS 10401 292 
VIMS 10404 237 
VIMS 10405 210 
VIMS 10406 195 
VIMS 10407 172 
VIMS 10408 114 
VIMS 10409 181 
USNM 1240011 230 


C. fulva 





CTAB (hexadecyltrimethylammonium bromide) was not 
added to the extraction and phenol was added immediately 
following incubation of the tissue at 37°C. DNA was pre- 
cipitated by the addition of 0.04 volume of 5M NaCl and 1.0 
volume of isopropanol. DNA was resuspended in 150 uL 
of sterile 0.1X TE (Tris-EDTA) and stored at —20°C. 


Morphological analyses 


Specimens of C. fulva and P. furcifer used in the mor- 
phological analysis were obtained from and measured 
at the National Museum of Natural History (Table 1). 
Ten putative hybrids were suitable for morphological 
analysis. The remaining five samples were not properly 
preserved and morphological analysis was not possible. 
Hybrid specimens as well as a small number of the C. 
fulva and P furcifer were frozen. Specimens obtained from 
the National Museum of Natural History were fixed in for- 
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malin. Morphometric and meristic characters were exam- 
ined as described in Smith (1971) by using dial calipers 
and a meter stick. The morphometric data were analyzed 
by using a sheared principal component analysis with a 
covariance matrix to confine the effect of size to the first 
principal component (Humphries et al., 1981; Bookstein 
et al., 1985; Stauffer et al., 1997). The meristic data were 
analyzed by using a principal component analysis with a 
correlation matrix. 


Allozyme analysis 


Horizontal starch ge] electrophoresis followed the pro- 
tecols described in Murphy et al. (1996). Gels (12%w/v; 
Starch Art Corp., Smithville, TX) were run on one of three 
buffer systems: Tris-citrate I] buffer (TC IF) (30 mAmps for 
14 hours), lithium hydroxide buffer (LIOH) (25 mAmps for 
14 hours), or Tris borate-EDTA buffer (EBT) (30 mAmps 
for 14 hours). Histochemical staining followed the proto- 
cols of Murphy et al. (1996), and locus nomenclature and 
allelic designations followed Shaklee et al. (1990). 

A preliminary survey of 16 loci in 15 individuals each of 
C. fulva and PF furcifer (Table 2) was performed to identify 
those loci for which the alleles were consistently differ- 
ent among the presumed parent species, C. fulva and P 
furcifer. All parent individuals and putative hybrids were 
then surveyed for all loci that demonstrated differences 
between the species. 


Nuclear DNA analysis 


An actin gene intron and the second intron in the S7 
ribosomal protein gene (Chow and Hazama, 1998) were 
investigated by using restriction fragment length poly- 
morphism (RFLP) analysis. Amplification primers and 
reaction conditions are listed in Table 3. The regions were 
amplified by using the PCR reagent system (GIBCO/BRL 
Life Technologies®, Bethesda, MD) and a 25-pL reaction 
cocktail (1X PCR buffer with MgCl,, 0.2 mM dNTP, 0.5 
uM primer, 2.5 U of Tag DNA polymerase, and 25-50 ng 
genomic DNA template). Some PCR reactions were per- 
formed with Platinum® Taq high fidelity (GIBCO/BRL 
Life Technologies®) with a 25-pL reaction cocktail (1X 
high fidelity buffer, 2 mM MgSO,, 0.2 mM dNTP, 0.2 pM 
primer, and 2.5 U of Platinum® Tag DNA polymerase high 
fidelity). In other cases, 1 uL dimethyl sulfoxide (DMSO, 
Fisher Scientific BP231-1, Pittsburgh, PA) was added to 
the reaction to increase sensitivity. 

The amplification products from two individuals of each 
species for both loci were digested with a panel of restriction 
enzymes to identify those that exhibited differences between 
the putative parent species. All samples were subsequently 
digested with those enzymes that demonstrated differences 
in the pilot study (Table 3). Digestion reactions (1.5-pL 10X 
buffer, 3 U restriction enzyme, and 4-uL PCR product) were 
incubated 2 to 18 hours at 37°C. Digestion products were 
separated on 2.5% agarose gels (1.25% Ultrapure Agarose, 
GIBCO/BRL Life Technology (R) + 1.25% NuSieve GTG (R) 
agarose, FMC Biochemical, Rockland, ME), stained with 
ethidium bromide and visualized under UV light. 
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Table 2 


Allozyme analysis: information includes loci, buffer sys- 
tems, and tissues. 





Locus Buffer Tissue 





Alcoho] dehydrogenase 
(ADH-1* 1.1.1.1) hver 


Creatine kinase (CK-B* 2.7.3.2) liver 
Creatine kinase (CA-C*2.7.3.2) liver 
Esterase, (EST-1* 3.1.1.1) liver 
Esterase, (EST-2* 3.1.1.1) liver 
Fumarase, (FH* 4.2.1.2) liver 


Glucosephosphate isomerase 
(GPI-A* 1.1.1.49) liver 


Isocitrate dehydrogenase 
UCDH-S* 1.1.1.42) liver 

Lactate dehydrogenase 
(EDITA 1.1.1.27) 


Lactate dehydrogenase 
(ELDER TILE liver 


muscle 


Malate dehydrogenase 
(MDH-A* 1.1.1.37) 


muscle 


Malate dehydrogenase 
(MDH-B* 1.1.1.37) 


Peptidase-B (PEP. 1” 3.4.11) liver 
Peptidase-S (PEP-2* 3.4.11) Hver 
Peptidase-C (PEP-3* 3.4.11) liver 


Xanthine dehvdrogenase 
(XDH* 1.1.1.204) liver 


muscle 


MtDNaA analysis 


The following regions of the mitochondrial genome were 
surveyed in the three putative parent species and hybrids; 
the adenosine 5’-triphosphatase subunit 6 (ATPase 6) 
gene, the 125/168 ribosomal RNA gene region, and the 
nicotinamide dehydrogenase subunit 4 (ND4) gene. Ampli- 
fication primers and reaction conditions are listed in Table 
3. Amplified products of two individuals of each putative 
parent species were screened with restriction enzymes to 
identify potential differences between C. fulva and P. fur- 
cifer (Table 3). All individuals were screened at the three 
regions with those enzymes that revealed differences in 
the preliminary study. Restriction digestion reactions were 
performed and visualized as described for nuclear DNA. 
For each individual, the haplotype designations of each 
region were combined in sequence creating a composite 
haplotype. 


Molecular data analysis 


Neis (1978) unbiased genetic distance was calculated 
from the allozyme data by using the computer program 
BIOSYS2 (Swofford and Selander, 1989). Mean nucleotide 
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Table 3 


PCR primers and conditions used in the RFLP analysis of mtDNA and nuclear intron regions. The forward primer is on top, the 
reverse is on the bottom. 





Region Primer sequence PCR conditions Citation Enzymes 


12SA-L: 5’ AAA CTG GGA TTA GAT ACC CCA CTA T 3° 
16SA-H: 5° ATG TTT TTG ATA AAC AGG CG 3’ 


125/165 94°C for 1 min., Palumbi et al., Ban Íl, Rsa I 
45°C for 1 min. 1991 


65°C for 3 min. 


95°C for 1 min., 
45°C for Imin., 
65°C for 3 min. 


95°C for 1 min., 
45°C for 1 min., 
65°C for 3 min. 


95°C for 30 sec., 
60°C for 1 min., 
72°C for 2 min. 
94°C for 1 min., 
48°C for 1 min., 
65°C for 3 min. 


H8969: 5’ GGG GNC GRA TRA ANA GEC T 3° 
L8331: 5’ TAA GCR NYA GCC TTT TAA G 3° 


ATPase 6 


Quattro! Dde I 


ARG-BL: 5’ CAA GAC CCT TGA TTT CGG CTC A 3’ 
LEU: 5’ CCA GAG TTT CAG GCT CCT AAG ACC A 3’ 


Bielawski and 
Gold, 1996 


BstO 1, Hpa UU, 
Mbo I, Rsa I 


S7 ribosomal S7RPEX2F: 5’ AGC GCC AAA ATA GTG AAG CC 3’ 
S7RPEX2R: 5’ GCC TTC AGG TCA GAG TTC AT 3’ 


Chow and 
Hazama, 1998 


Alu I, Dra I 


F3: 5’ ATG CCT CTG GTC GTA CCA CTG G 3’ 
R1: 5’ CAG GTC CTTACG GAT GTC G 3° 


Actin intron Cordes, 2000 Hinfl 


? Quattro, J 1999 Personal commun. Department of Biology, Univ. South Carolina, Columbia, SC 29208. 





sequence divergence was calculated for the mito- 
chondrial DNA RFLP results by using the equation 
of Nei and Li (1979) for fragment data with weight- 
ing based on Nei and Tajima (1983) as performed by 
the computer program REAP (McElroy et al., 1992). 
Genetic distance was not calculated for the nuclear 
DNA loci owing to an absence of shared restriction 
fragments between C. fulva and P. furcifer. 


e C. fulva 
m P. furcifer 
Putative Hybrid 


Results 


-0.6 -0.4 -0.2 0 02 0.4 0.6 


Morphological analysis Sheared PC 2 


Of 44 counts and measurements analyzed, mean 
values of 38 were different between C. fulva and P. 
furcifer, and the putative hybrids were intermediate 
in morphology to presumed parents for 33 of these 
characters (Table 4). The mean value for the puta- 
tive hybrids exceeded the values of either presumed 
parent species for four characters (caudal peduncle 
scales, interorbital width, anal base length, and 
pectoral length) and was lower than either pre- 
sumed parent species for two characters (orbit length and 
depressed dorsal length). In the principal component analy- 


Figure 2 


Graph of principal component 2 for the morphological data and 
factor 1 for the meristic data for ten Cephalopholis fulva, ten 
Paranthias furcifer, and ten putative hybrids. The arrow indicates 
putative hybrid H7, an individual identified as a backcross to C. 
fulva, according to the genetic data. 





component from the morphological analysis (variation in 
shape) and the first factor of the correlation matrix of the 


sis, the first component accounted for 66% of the variation 
in the data and the second accounted for 23%. Suborbital 
width, caudal peduncle length, caudal peduncle to upper fin 
rays and caudal penduncle to lower fin rays were the char- 
acters with the highest loadings. In the meristic data, the 
first principal component accounted for 62% of the varia- 
tion. Dorsal rays, gill rakers, and transverse scale rows 
had the highest loadings. A plot of the second principal 


meristic data shows a discrete difference between the puta- 
tive parents (Fig. 2). Hybrid individuals were shown to be 
intermediate between the putative parent species. 


Allozyme analysis 


Sixteen allozyme loci were surveyed in 15 C. fulva and 15 
P. furcifer to identify those that had different alleles in 
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Table 4 
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Ranges of counts and measurements for ten Paranthias furcifer, ten Cephalopholis fulva, and ten putative hybrids. Measurements 
are given in millimeters. Raw measurements were divided by the standard length and multiplied by 1000. Means are given in 
parentheses, L.L. = lateral line, and the measurement for one hybrid discarded because of a broken third dorsal spine (*). 





Measurement 





C. fulva 





dorsal rays 

anal soft rays 

pectoral rays 

gill rakers 

scales above the L.L. 

scales below the L.L. 
transverse scale rows 
caudal peduncle scales 
head length 

head width 

head depth 

snout length 

suborbital width 
interorbital width 

orbit length (diameter) 
postorhital head length 
maxillary length 

lower jaw length 

snout to angle of preopercle 
maxillary width 

tip of lower jaw to gular notch 
body width 

body depth 

caudal peduncle depth 

tip of snout to dorsal origin 
tip of snout to pectoral base 
tip of lower jaw to pelvic base 
dorsal base length 
depressed dorsal length 
anal base length 

depressed anal length 

end of dorsal to caudal hase 
length of caudal peduncle 
pectoral length 

pelvic length 

dorsal spine I length 
dorsal spine ITI length 
dorsal spine IX length 

anal spine I length 

anal spine IT length 

anal spine 111 length 
caudal base to tip of upper rays 


caudal base to tip of middle rays 


caudal base to tip of lower rays 


TX, 15-16 (16) 


8-10 (9) 
32-35 (34) 
22-28 (25) 

6-8 (8) 
22-27 (25) 
64-84 (71) 
40-49 (46) 

374-427 (402) 
173-227 (200) 
248-288 (266) 
86-122 (107) 
44-52 (47) 

64-78 (71) 

62-76 (68) 

217-2483 (232) 
166-194 (182) 
176—200 (186) 
257-306 (282) 
42-56 (48) 

119-181 (145) 
154-223 (178) 
325-385 (357) 
128-139 (134) 
386—422 (407) 
287-415 (371) 
398-449 (420) 
523-549 (538) 
609-663 (629) 
166-187 (174) 
268-320 (294) 
129-150 (141) 
166-198 (178) 
257-292 (274) 
186-211 (197) 
54-66 (62) 

103-138 (123) 
97-150 (127) 
52-72 (62) 

95-121 (106) 
78-118 (108) 
197-246 (216) 
204-248 (222) 
197-242 (217) 





Putative hybrids 


P. furcifer 





IX, 17-18 (17) 
9-10 (9) 
35-37 (36) 
24-35 (31) 
10-12 (11) 
26-30 (28) 
75-91 (84) 
44-52 (47) 
320-353 (335) 
150-189 (164) 
212-250 (233) 
79-111 (94) 
29-35 (32) 
73-89 (82) 
56-67 (62) 
173-210 (192) 
120-163 (141) 
117-149 (138) 
213-246 (229) 
30-45 (41) 
90-174 (123) 
155-181 (165) 
302-361 (332) 
118-145 (129) 
334-373 (354) 
293-340 (319) 
359-420 (382) 
532-561 (543) 
574-648 (609) 
171-192 (182) 
248-273 (261) 
148-165 (157) 
194-217 (203) 
253-295 (279) 
143-192 (180) 
43-70 (61) 


90-119 (106) ` 


104-123 (111) 
35-66 (53) 

91-112 (101) 
97-121 (106) 
228-304 (280) 
164-192 (175) 
278-318 (294) 


IX, 18-19 (19) 


9-11 (10) 
36-39 (38) 
32-39 (36) 
10-14 (12) 
26-32 (30) 
85-96 (91) 
43—48 (46) 

259-301 (280) 
130-162 (142) 
183-241 (212) 
53-79 (67) 
20-24 (22) 
76-89 (80) 
52-89 (65) 
146-171 (159) 
96-115 (104) 
98-114 (105) 
176-200 (188) 
28-38 (34) 
69-121 (89) 
138-162 (150) 
282-354 (316) 
102-122 (111) 
321-361 (329) 
265-293 (277) 
321-384 (351) 
543-613 (583) 
603-658 (635) 
165-204 (180) 
232-278 (253) 
145-175 (162) 
174-241 (219) 
249-290 (274) 
150-187 (173) 
40-64 (53) 
86-117 (105) 
76-105 (9-4) 
36-48 (43) 
78-98 (86) 
74-101 (87) 
321-373 (347) 
122-143 (128) 
308-346 (322) 





the two species (Table 2). At four loci (CK-B*, FH*, LDH- 
B*, and ICDH-S*) C. fulva and P. furcifer had different 
alleles. Forty C. fulva, 28 P. furcifer, one C. cruentata, and 
ten putative hybrids were subsequently screened at the 
four diagnostic allozyme loci (Table 5). Eight of ten puta- 


tive hybrids were heterozygous at all four diagnostic loci. 
One putative hybrid was heterozygous at all loci except 
the LDH-B* locus, for which it displayed two alleles 
characteristic of P. furcifer, and another individual was 
heterozygous at all diagnostic loci, except the FH* locus 
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Table 5 


Allozyme genotypes of Cephalopholis fulva, Paranthias furcifer, C. cruentata, and putative hybrid individuals at four diagnostic 
loci, CK-B*, FH*, ICDH-S*, and LDH-B*. n = number of fish in a sample. 


CK-B* 


C. fulva *100/100 
P. furcifer : *50/50 
*100/100 


*50/100 
*50/100 
*50/100 


C. crucntata 


Putative hybrids 


Table 6 


RFLP genotypes for the nuclear intron and mtDNA loci. 
Some bands (*) were inferred so that fragments would sum 
to the total (uncut) size of the amplified gene region. 


Approximate 

Locus Enzyme Allele size (bp) 

400, 50° 
450 


Sip, 50075 
1200 
575.525. 151, 25. 
600, 50° 
360, 150, 90, 50° 
650 


600. 500, 300, 200 
600, 300, 250, 250, 200 
450, 375, 300, 
275, 100, 100° 


1600 
1100, 500 


1250, 625, 25° 
1250, 450, 200 
1450, 450 


1400, 500 
1900 
1000, 900 


700, 500, 400, 300 
550, 500, 400, 300, 150 
525, 400, 300. 275, 
250, 80, 70 


1025. 500, 375 
610, 500, 400, 390 
500, 400, 350, 
300, 300, 50° 
900, 425, 325, 250 


actin intron Hinf I 


S7 intron Dra } 


ATPase 6 Dde } 


128/168 


QW Or BOP we 


Ban 1} 


BstO H 


Hpa H 


Mbo } 


awr Jarr Or OY 





for which it displayed two alleles diagnostic of C. fulva. 
Cephalopholis cruentata was distinguished from C. fulva 
and F. furcifer at the ICDH-S* locus and was eliminated as 
a potential parent species. 





FH* ICDIL-S' 


*90/90 *95/95 *100/100 
*100/100 100/100 “75/75 
*90/90 105/105 75/75 


*90/100 *100/95 *75/100 
*90/100 100/95 *75/15 
*90/90 "100/95 *75/100 


LDII-B 





Table 7 


Genotypes of Cephalopholis fulva, Paranthias furcifer, C. 
cruentata, and putative hybrids for the short actin intron 
and the second intron in the S7 ribosomal protein region. 
n = number of fish in a sample. 


Actin intron S7 intron 


Species Hinf1 Dra } 


Cephalopholis fulva f AIA 


Paranthias furcifer 
C. cruentata 


Putative hybrid 





Nuclear intron regions 


An actin intron approximately 450 base pairs in length 
was amplified. The region was surveyed with 14 restric- 
tion enzymes, of which Hinf I showed a genetic difference 
between species. After digestion with Hinf I, all C. fulva 
demonstrated allele A, and all P furcifer diplayed aliele B 
(Tables 6 and 7). All fifteen putative hybrids were hetero- 
zygous for both alleles. 

The second intron region of the S7 ribosomal protein, 
which was approximately 1200 base pairs in length, was 
screened with thirty-five enzymes. Two enzymes, Dra I 
and Alu I, demonstrated differences between P. furcifer 
and C. fulva. Paranthias furcifer and C. cruentata both 
exhibited allele D after digestion with Dra I (Tables 6 and 
7), Cephalopholis fulva was variable at this locus—fifty in- 
dividuals were homozygous for allele A, three homozygous 
for allele E, and four heterozygous for alleles A and E. All 
fifteen putative hybrids were heterozygous at this locus 
with one of the C. fulva alleles (A or E), and the P. furcifer 
allele (D). Digestion of the second intron in the S7 region 
by Alu I produced a large number of small fragments that 
were not easily interpreted and the data were not used 
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to identify hybridization between C. fulva and P. furcifer. 
With this enzyme, however, C. cruentata had a unique al- 
lele and was thus eliminated as a putative parent for the 
hybrid individuals. 


Mitochondrial DNA 


Allelic differences between C. fulva and P. furcifer were 
found in all three mitochondrial gene regions. The ATPase 
6 region was screened with six enzymes, one of which, Dde 
I, showed differences between C. fulva, P. furcifer, and 
C. cruentata (Tables 3 and 6). The 125/168 region was 
screened with seven enzymes. Two of these, Ban I] and Rsa 
I, demonstrated differences between C. fulva, P. furcifer, 
and C. cruentata. The ND4 region was screened with nine 
restriction enzymes, four of which (BstO I, Hpa H, Mbo I, 
and Rsa 1) showed differences between the species; mtDNA 
composite haplotypes were unique to each species (Table 
8). All 15 putative hybrids in the study had a composite 
haplotype matching the common haplotype of C. fulva, 
indicating that it was the maternal parent for all hybrid 
individuals. Three Epinephelus guttatus specimens were 
screened at the ND4 and ATPase 6 regions and showed a 
unique composite haplotype; therefore this species was not 
included in the study as a putative parent species. 


Discussion 


Morphological analysis 


In most cases, F, hybrids should be morphologically 
intermediate to the parent species and have low variation 
within characters among themselves. Backcross individu- 
als, because of random sorting of chromosomes, should 
have higher variation within intermediate characters and 
could fall anywhere in the morphological range of the pure 
parent species (Anderson, 19-49). In a principal component 
analysis plot, a backcross hybrid’s score would be expected 
to be closer to the parent species to which the hybrid 


Table 8 
Composite haplotypes of Cephalopholis fulva, Paranthias 
furcifer, and C. cruentata for the mitochondria] DNA data. 
Haplotypes are given for the following sequences of mt- 
DNA loci and enzymes: 1) ATPase 6—Dde I; 2) ND4—BstO 
I, Hpa Il, Mbo 1, Rsa I; and 3) 128/165—Rsa I, Ban Il. n = 
number of fish in a sample. 





Species 


Composite haplotype 





AAABBAA 
AAABCAA 


BBCADBB 
BBDADBB 
CCDCACB 


AAABBAA 


Cephalopholis fulva 


Paranthias furcifer 


C. cruentata 


Putative hybrids 





backcrossed, whereas an F} hybrid’s characters would be 
expected to be in the center, closer to an average of the 
scores of the parent species. 

A plot of the second principal component of the morpho- 
logical analysis and the first factor of the meristic analysis 
(Fig. 2) shows that C. fulva and P. furcifer are well segre- 
gated according to morphological characters. The putative 
hybrid individuals were clustered in between the parent 
species. The post-F, hybrid detected by using genetic anal- 
yses, indicated by an arrow in Figure 2, clustered with the 
putative F, hybrid individuals. 


Genetic analyses 


Results of the allozyme analysis also supported hybridiza- 
tion between C. fulva and P. furcifer in Bermuda. The puta- 
tive hybrids were heterozygous at four distinguishing loci, 
with the exception of two individuals that were homozygous 
at one diagnostic locus each. One individual was homozy- 
gous at the LDH-B* locus for the *75 allele. All twenty-eight 
P. furcifer were homozygous for this same allele, indicating 
hybrid backcrossing to P furcifer. Another hybrid individual 
was homozygous at the FH* locus for the *90 allele. The 
40 C. fulva sampled in this study were homozygous for 
this allele, suggesting hybrid backcrossing to C. fulva. This 
hybrid individual was among those included in the morpho- 
logical study and, as shown in the principal component plot 
(Fig. 2), was morphologically the most similar to C. fulva. 
It was not possible to distinguish F, hybrids and backcross 
hybrid individuals and henceforth the two individuals 
described above are referred to as post-F, hybrids. Because 
all members of the presumed parent samples were homozy- 
gous at all loci for diagnostic alleles, there was no evidence 
of introgression between C. fulva and P furcifer. 

The nuclear intron data were consistent with the allo- 
zyme data and supported the hypothesis of hybridization 
between C. fulva and P. furcifer. Because all hybrid indi- 
viduals were heterozygous, post-F, hybridization was not 
evident at these loci. Alleles present at both nuclear DNA 
loci were unique between parent species and there was no 
indication of introgression between these species. 

The mtDNA data clearly showed that C. fulva was the 
maternal parent for all putative hybrids, including the two 
post-F, hybrids. This finding suggests a strong gender bias 
in hybridization. All C. fulva had composite haplotypes quite 
distinct from those of P furcifer, and there was no evidence 
of mtDNA introgression between the two parent species. 

Overall, the genetic and morphological analyses sug- 
gest that all but two of the 15 putative hybrids were F} 
individuals representing first generation hybridization 
between a female C. fu/va and a male P. furcifer (Table 9). 
The occurrence of two post-F, hybrids indicates that F, 
hybrids are fertile, and the genotypes of the two post-F, 
hybrids demonstrate that F, hybrids can backcross with 
either parent species. 


Hybridization in Bermuda 


Hybridization between C. fulva and P. furcifer is known 
from only certain localities in the tropical Atlantic, despite 
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Table 9 


Classification of all putative hybrid individuals as F} and post-F, individuals based on morphological and genetic data. 





Sample ID Classification 


H1 post-F, 
H2 

H3 

H4 

H5 

116 

H7 

H8 

H9 

H10 

H11 

H12 

H13 

H14 

H15 

USNM 124001 


gai iezh ieo ea ee al E iel 





broad overlap in the geographical ranges of the two spe- 
cies. Bermuda exists at the northern range of C. fulva and 
P. furcifer, and the two species have restricted spawning 
times. Cephalopholis fulva and P. furcifer both spawn in 
Bermuda from May to early August (Smith, 1958, as cited 
in Thompson and Munro, 1978), and spawning individuals 
of the two species have been sampled in the same location 
(Burnett-Herkes, 1975; B. Luckhurst, personal observ.). 
Similarly, two hybrids have been reported from Jamaica 
where the parent species also have overlapping spawning 
times (Thompson and Munro, 1978). 

Hybridization between C. fulva and P. furcifer could 
be the result of directed interspecific interactions, or the 
chance meeting of gametes spawned at the same time 
in the same general location. However, because all 13 F} 
hybrids were the result of C. fulva eggs fertilized with P. 
furcifer sperm, it appears that there is a gender bias in 
hybridization. Possible hybridization scenarios include 
differences in sex ratio between C. fulva and P. furcifer, 
a biochemical block on fertilization of P furcifer eggs by 
C. fulva sperm, or “sneaker” P. furcifer males in C. fulva 
spawning groups. 

Based on the number of individuals collected in Ber- 
muda over the last two years, hybridization between C. 
fulva and P. furcifer seems to be relatively rare. Although 
the reproductive status of the hybrids is unknown, the oc- 
currence of two fertile female F, hybrids (one of which was 
ripe), a spent male F, hybrid, and the presence of post-F, 
hybrids, it can be concluded that some F} hybrids are ca- 
pable of producing viable offspring. However, reproduction 
of F, hybrids does not appear to be very extensive. Within 
the limits of the samples analyzed there was no evidence 
of introgression of alleles between parent species. 


Data 


allozyme, mtDNA, nDNA 

mtDNA, nDNA 

mtDNA, nDNA 

mtDNA, nDNA 

mtDNA, nDNA 

mtDNA, nDNA 

morphology, allozyme, mtDNA, nDNA 
morphology, allozyme, mtDNA, nDNA 
morphology, allozyme, mtDNA, nDNA 
morphology, allozyme, mtDNA, nDNA 
morphology, allozyme, mtDNA, nDNA 
morphology, allozyme, mtDNA, nDNA 
morphology, allozyme, mtDNA, nDNA 
morphology, allozyme, mtDNA, nDNA 
morphology, allozyme, mt DNA, nDNA 
morphology 


Intergeneric hybridization 


Intergeneric hybridization in animals seems to be rela- 
tively rare. It is believed that the ability to hybridize is an 
indication of evolutionary relatedness and that divergent 
taxa should have lost the ability to interbreed through the 
evolution of reproductive isolating barriers (Sibley, 1957). 
However, there are several examples of intergeneric hybrid- 
ization in fishes, most notably in the cyprinids (Hubbs, 1955; 
Smith, 1973; Aspinwall et al., 1993; Stauffer et al., 1997). 

Intergeneric hybridization between two such ecologi- 
cally different species as C. fulva and P furcifer has also 
been noted in the lutjanids. Poey (1860) described an in- 
tergeneric hybrid between Lutjanus synagris and Ocyurus 
chrysurus. Both Loftus (1992) and Domeier and Clarke 
(1992) presented reviews of this case of hybridization and 
Loftus (1992) theorized that the species were capable of 
interbreeding on the basis of overlap of spawning time 
and habitat. Both Loftus (1992) and Domeier and Clarke 
(1992) surmised that O. chrysurus should not constitute 
its own genus but be included within the genus Lutjanus, 
owing to its ability to hybridize with a member of that 
genus. Further support for this revision was provided by 
Chow and Walsh (1992), who reported a high genetic simi- 
larity between O. chrysurus and several Lutjanus species 
of the western North Atlantic. They suggested that the 
species appeared to be misplaced due to its unique mor- 
phological features. 

A parallel situation exists with the epinepheline serrands 
because the genetic distance data suggest that P. furcifer 
does not belong in a separate genus from Cephalopholis. 
Recently, Craig et al. (2001) used mitochondrial DNA se- 
quencing to demonstrate a close phylogenetic relationship 
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between the genera Paranthias and Cephalopholis. They al- 
so cited several morphological and ontogenetic similarities 
between the two genera. In the present study, an analysis 
of 16 allozyme loci produced a Nei’s (1978) unbiased genetic 
distance of 0.356 between P furcifer and C. fulva. which is 
within the range reported between Ocyurus chrysurus and 
species of Lutjanus (Chow and Walsh, 1992). The mitochon- 
drial DNA sequence divergence data indicated that P fur- 
cifer is closer to, or at least no farther from, C. fulva (0.036) 
than is C. cruentata (0.042) and that the congeners are most 
distant, (0.052). Both values are within the range reported 
between other congeneric serranids (Graves et al., 1990). 

Taxonomic revision of P. furcifer is not recommended at 
this time; however, the incidence of hybridization and the 
small genetic distance between C. fulva and P furcifer sug- 
gest that a full-scale phylogenetic analysis of the subfam- 
ilv is warranted. Such a study would have to include more 
species of Cephalopholis, other members of the subfamily 
Epinephelinae, and P colonus, the eastern Pacific gemi- 
nate species of P. furcifer. 
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